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We use a semiclassical approach for analysing the tunneling transport through 
a normal conductor in contact with superconducting mirrors. Our analysis of the 
electron-hole propagation along semiclassical trajectories shows that resonant trans- 
mission through Andreev levels is possible resulting in an excess, low-energy quasi- 
particle contribution to the conductance. The excess conductance oscillates with 
the phase difference between the superconductors having maxima at odd multiples 
of 7r for temperatures much below the Thouless temperature. 
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1. Introduction 

Recent experimental and theoretical work on diffusive charge transport in mesoscopic N/S samples 
[1-39] have revealed a strong energy dependence of an excess quasiparticle contribution to the low 
temperature conductance of normal (N) parts in close proximity to superconductors (S) (for a review 
see e.g. A characteristic energy scale is set by the Thouless energy Etk = TiD/ L 2 S below which 

a re-entrance to normal conduction is seen as the bias voltage or temperature is lowered (D is the 
diffusion coefficient, L$ is the distance between the superconductors). In samples of the Andreev 
interferometer type with two N/S interfaces the excess conductance oscillates as a function of the 
phase difference <j> between two superconductors. Conductance maxima occur at even multiples of 
7r, their magnitude peaks at Etk and becomes vanishingly small at low energies. These oscillations 
have been explained by Nazarov and Stoof as a "thermal effect" [M |j| and can be understood in 
terms of competing contributions from the condensate wave functions to the density of states or in 
terms of quasi-particle trajectories between the N/S interfaces 

In Ref. ||, a strong interference effect due to resonant transmission of quasiparticles through 
Andreev levels was shown to take place in superconductor-normal metal-superconductor (S/N/S) 
heterostructures at temperatures corresponding to energies much below the ballistic Thouless energy 
Etk = ^vf/Ls (vp is the Fermi velocity). The effect results in giant peaks in the conductance — 
proportional to the number Nj_ of conducting transverse modes — whenever the phase difference 
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4> between the superconductors is an odd multiple of n. Arguments were presented in Ref. |Hj for 
this result to be valid also in the diffusive diffusive transport regime. 

Andreev levels are bound states formed in a normal sample clement by successive Andrccv reflec- 
tions of quasiparticle excitations at two S/N- interfaces. A peculiar feature of these bound states is 
that they carry the supercurrent (if any) between the two superconductors. These levels can also 
carry a normal transport current if the sample is coupled to reservoirs of normal electrons p!o| . Such 
a transmission is of a resonant type if the coupling is weak enough not to destroy the Andreev levels 
themselves. This is the case when the reservoirs are coupled through tunneling barriers of low trans- 
parency. It is important that, in addition, such barriers serve as quantum scatterers extended in two 
dimensions; they split quasiclassical electron trajectories incident from the reservoir, or returning 
towards the reservoir after having been Andreev reflected from an N/S boundary. This enables the 
trajectory of a quasi-particle with zero excitation energy (measured from the Fermi level) to connect 
both S/N- interfaces (necessary for picking up information about cf> = cf>i — cf>2, 4>\ and <j>2 being the 
phases of superconductors 1 and 2 ) and a reservoir (to affect the current) as shown in Fig. 1. 

In this work, which is based on an analysis of the quasiparticle trajectories in a disordered normal 
conductor weakly coupled to reservoirs, we calculate how the two-terminal conductance between 
two normal resrvoirs depend on both the superconductor phase difference and temperature. We 
pay special attention to non- Andreev (normal) reflections at the N/S boundaries. Including these, 
the character of the resonance is changed shifting the maximum of the amplitude of conductance 
oscillations to nonzero temperatures. In our concluding section we discuss recent experimental results 
of Ref. EL where an anomalous low-temperature behavior of the conductance was observed. 



Figure 1 shows the geometry of the sample under consideration. A mesoscopic normal region in 
contact with two superconductors is coupled to two reservoirs of normal electrons through tunneling 
potential barriers of low transparency. The quasiparticle motion in the normal region is assumed to 
be affected by a smoothly varying disorder potential and is treated semiclassically. 

According to the Landauer-Lambert formula || the conductance of the system corresponding to 
a current between two normal electron reservoirs can be written in terms of the probability R a (E) 
for electrons impinging from one reservoir with energy E to be reflected as holes back into the same 
reservoir and in terms of the probability T (E) for them to be transmitted to the other reservoir as 
electrons. For the case of two equal barriers one has 



where fo(E) is the Fermi distribution function. Our objective is to calculate the "electron- hole 
transmission" probability R a {E) in Eq. (p. This is sufficient, since it can be shown that due to the 
destructive interference among the relevant trajectories, Tg(E) is only weakly affected by the super- 
conductor phase difference </>. Hence the effective transmission probability T{E) for our problem, to 
be calculated below, is given by R a {E). 

An electron impinging on the lower potential barrier in Fig. 1 from the reservoir gets, after possibly 
interacting with the barrier and the two S /N interfaces in the N-region, backscattered in the electron 

and hole channels with probability amplitudes Cp 6 ^(r||) and Cp^(r||) respectively. The coordinate 
rii = (0, y, z) lies in the plane of the barrier, the x-axis being perpendicular to the barrier plane. 
Accordingly, the wave function of the electron u(r||) and hole i>(i"||) in the reservoir can be written 



2. Formulation of the problem 




(1) 



Superlattices and Micro structures, Vol. ??, No. ?, 1999 



3 



RESERVOIR 2 




Fig. 1. (a) Structure analyzed in the text, with superconducting (S) and disordered normal (N) elements labeled. 
The normal element is coupled to two reservoirs through tunnel barriers shown as thick black lines. Semiclassical 
trajectories of electron- and hole-like quasiparticle excitations are indicated by solid and dashed lines, respectively. 
Sections of the trajectories going between N/S boundaries 1 and 2 are connected by Andreev and normal reflections 
at the points shown as black dots. An electron incident from the resorvoir (solid up-arrow) may be reflected back 
either as an electron (solid down-arrow) or as a hole (dashed arrow), (b) Sequence of scattering events along the 

trajectory shown in (a). L n is the length of the trajectory in section n; a^'^ are the amplitudes of the corresponding 
electron and hole excitations at the Fermi energy (n = 0, ±1, . . .). 



as 
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The probability for electron-hole transmission, which we need to calculate, is connected with the 

II 
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amplitude Cp (rii) by the relation 
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(3) 



where S is the area of the tunnel barrier. 

The wave function in the normal region N near the potential barrier is also characterized by 
electron and hole components, 



u(r) 
v{r) 



= r*P|l r ll 



4 e, (p|)e*i-'» + C('|)e' 



Bi h) (r {] )e- ip * )x 



(4) 



A slowly varying potential in the N region (on the scale of the Fermi wavelength Xp) is responsible 
for the semiclassical nature of the wavefunction (||) and the slow spatial variation of the factors 
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A( e > h '(ru) and B^ e < h '(rn). The matching conditions at the barrier, which can be expressed in terms 
of a unitary 2x2 matrix, describe the coupling between C^ e ' h \r\\) and A^ e ''^(r||), B^ e ' h ^ (r|| ). We 
assume that these functions smoothly — in the semiclassical sense — go to zero at the perimeter 
of the injector reservoir and are equal to zero in the plane x = outside the injector. 

Using the language of semiclassical propagation we now construct the wavcfunction in the N- 
region (see Fig. 1) for the electron impinging on the lower barrier from the reservoir. We do so by 
mapping it to the wave function (^) at the injector barrier. Such a mapping is possible if the disorder 
potential does not cause any noticable divergence of a tube of trajectories with the transverse size 
of the order of as they propagate for a certain length. In our case this characteristic length is 
the length covered by quasiparticle diffusing across the sample (see, e.g., ]42[|). 

The procedure for doing the semiclassical mapping can be reduced to the following. For a given 
point r in the N region we have to introduce all different trajectories which connect this point 
with points on the barrier (where the momentum of the quasiparticle is p), and experiencing all 
possible sequences of scatterings induced by N/S boundaries and barriers. According to Ref. |42] |. the 
semiclassical wavefunction can be constructed as a sum of partial contributions corresponding 
to these trajectories, expressed in terms of the classical action Si = J L pdl. As a result one has 

u(r) = J2 A i e) ^ elS ' C>/h 
i 

w(r) = j2 A i h)( - r y s ' h>/h - ( 5 ) 

i 

If a trajectory is split by interacting with a tunnel barrier or N/S interface, the wavefunction along 
that trajectory undergoes a transformation described by the scattering matrix already mentioned. 
The smoothly varying function A, (r) can be found from the continuity equation for the current 
density, which together with the Hamilton- Jacobi equation for the action, guarantees that Eq. (^) is a 
solution of the Schrodinger equation |f42f . Furthermore, one can readily verify that the wavefunction 
, constructed as a sum over trajectories, satisfies the boundary conditions corresponding to the 
scattering at the barriers and interfaces. 



3. Electron-hole transmission at low energies 

The wave function formally constructed in Eq. (^|) does not permit us to carry out concrete cal- 
culations in a general case. The reason lies in the complications that arise due the bifurcations of 
trajectories as they undergo Andreev and normal reflections at the N/S boundaries. The situation 
is drastically simplified in the region of small energies, E <C Erh, where reflections in the Andreev 
channel send the quasiparticle back along the incident trajectory. This implies that the trajectory bi- 
furcations disappear and the problem reduces to an analysis of one-dimensional quasiparticle motion 
along a single trajectory with centers for back (Andreev) and forward (normal) scattering. 

A trajectory of this type is shown in Fig. 1. In this case the problem to find the wavefunction 
at the boundary [see Eq. (J2J) ] is reduced to a quantum scattering problem for the configuration 
shown in Fig. lb. The points of reflection at the N/S boundaries and the points of scattering at 
the tunnel barrier are shown with black dots and black bars, respectively. Propagation between 
these points is coherent in both electron and hole channels, which is illustrated by dashed and solid 
lines of equal lengths. For quasiparticles with finite energy E the phase gains along the electron 
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and hole trajectories do not completely cancel since the momenta are now different. The resulting 
decompensation effect is of order E/Exh- f\ 

Wave functions in the adjacent sections of Fig. lb are connected by a scattering matrix describing 
Andreev and normal reflections at the N-S boundary. The problem of evaluating the electron-hole 
transmission can now be formulated as follows. An electron in the reservoir arrives at the tunnel 
barrier (solid up-arrow in Fig. lb), and we have to find the amplitude of the outgoing wave function in 
the hole channel (dashed up- arrow). This problem is reduced to solving a set of matching equations 
for amplitudes of electron an and hole excitations in every section of propagation between 
scattering points. A further simplification follows as a consequence of the resonant transmission 
caused by multiple electron-hole transformations. Such a resonance occurs, as we will show, if the 
superconductor phase difference 4> is close to an odd multiple of it corresponding to a large number 
of trajectories contributing to a constructive interference between scattering events. If |rjv| <C 1, 
where r^ is the probability amplitude for non- Andreev (normal) reflection at an N/S boundary, the 
main contribution comes from trajectories wich do not include successive reflections at the same N/S 
boundary. Taking these observations into account the following set of algebraic matching equations 
emerge, 
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The amount of phase gained after propagation along the trajectories in section n is 

(E) = jf P ^ (E)di/h « $ n ± r„£;/?i , 

where $„ = < &^ e, ' l ' ) (0) and r n is propagation time in section n. The quantities r^' 2 ^ and r^' 2 ^ are, 
respectively, the probability amplitudes for Andreev and normal reflection at N/S boundaries 1 and 
2. Phases and amplitudes of an electron or hole along the semiclassical paths arc defined in such 
a way that no phase has been gained at the beginning of a particular electron or hole section n. 
Hence the amplitude is a^' h ^ at the beginning of the section, and a^'^e^*™ h) at its end. The 
phase gain between the tunnel barrier and the left N/S boundary (N/S boundary 1) is denoted by 
^q 6 ' 1 -*. In Eq.(^) a coefficient e r <C 1 characterizing the coupling through the tunnel barrier has 
furthermore been introduced. We note that one can show that the large phases $ n can be removed 
from the set of equations (^). This is a manifestation of the fact that the electron and hole phase 
gains compensate each other at E = 0. 

According to our construction the probability |C p (r||)| 2 of an electron-hole transmission at the 

point r || of the barrier is related to the amplitude on the trajectory of Fig. 1 corresponding to 



f The semiclassical trajectories for an electron of energy E and a reflected hole of energy —E are to be 
considered identical since they separate by less than Af while diffusing a distance Ls if E -C Eth- 
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injection at point r 1 1 as, 

|(7 p (r||)| 2 = e r |aW| 2 . (7) 

Equations (^) and (Q) together with give the complete solution for the oscillatory, (^-dependent 
part of the excess conductance. 

When e r = 0, Eq. (|J) gives a set of Andreev levels E a . Therefore, one can expect that if e r <C 1 
the trasmission probabilty | Cp /i '' > (r || ) | 2 is of the Breit-Wigner form. Indeed, as shown in Appedix 1, 
the transmission probability |Cp /l ' ) (r || ) | 2 in this limit can be expressed as 

|C P w ( r|l )l ~ A{E-E a y + elEf (8) 

Here E is the spacing of Andreev levels generated by the trajectory when e r = and n = 0, a^ a 
is the solution of Eq. (||) for E = E a and e r — 0. Because of the random variation of propagation 
times t„, the functions Oq 2 are localized along the one-dimensional ladder shown in Fig. lb. 

Further simplifications arise as one integrates the electron-hole transmission probability | C p (r || ) | 2 
[see Eq. (||)] over the area of the injector area and over the the Fermi surface in momentum space. 
This integration corresponds to summing over different trajectories and one can think of it as aver- 
aging over various distributions of r„. It follows that T(E) can be expressed as 




where N± — Spj^/h, S is the area of the injector and ((■■•)) implies an averaging over r n . If 
ksT 3> e r E() one can neglect the width of the resonance for relevant energies, E ~ kT, and conclude 
that 

T{E) = N^(^e r \a { ^8{[E E a ]/E )\\ . (10) 

The distribution of propagation times r„ depends on the details of the disordered potential in 
the mesoscopic normal region. These are not known, but it is natural to assume that propagation 
times along different sections of the semiclassical trajectory (see Fig. la) are uncorrelated. Under 
this assumption one can, as detailed in Appendix 1, directly express the transmission probability 
T(E) in terms of the average density of Andreev states coupled with the reservoir jj] 

T(E) = N x E Th ((v(E))) (11) 

with 

E Th = {{E )) = h((r- 1 )) . (12) 

Here v(E) is the density of Andreev states generated by a given semiclassical trajectory when e r = 0. 

In order to proceed with an analytical approach we choose a Lorentz form for the distribution 
function P(r), 

P ( r ) = -7 4f^' ( 13 ) 

f It is necessary to use the fact that ( ( \a^a \ 2 5(E — E a ) \ \ does not depend on n if the r n 's are uncorrelated 
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As shown in Appendix 2, this choice permits us to derive an analytical expression for the averaged 
density of state 

oo 
— oo 

(we have used 7 = r, see JhJ). The quantity vq(E) in the integrand of eq. ( |l~4| ) is the density of 
states of the periodic chain of Fig. lb, if r„ = f for all n. It is straightforward to find this density 
of states to be 

Vo (E) = - — = ^ , (15) 

The energies E min and E max are lower and upper edges of the energy band of the periodic chain of 
Fig. lb. They can be expressed in terms of the "hopping integral" r 1 ^ as 

£ min>roax = E Th Js<j> 2 + (\r$\ ± |r^|) 2 <50 = Min|0-7r(2fc + l)|. (16) 

Finally, using Eqs. (0), ([ll] and[IJ) one finds that the low-temperature conductance can be expressed 
as 

2e 2 e r V2 
h T~ Jo cosh 2 (a;/2T) 

{ \/(4:e 4 4- E min )(4:X 4: 4- ^tm) 4- E min E max 
(4* 4 + E^ in )(4x* + E^ ax ) 

where T = T/Tth with T T /» = E Th /k B defined by Eq.@. 

The excess conductance G(0) is plotted as function of phase difference in Fig. 2, while the maxi- 
mum oscillation amplitude G max (T) is plotted as a function of temperature in Fig. 3. We emphasize 
two distinguishing features: 1) the excess conductance has sharp maxima at <f> — n(2k + 1), k = 
0,±1, 2) The peak oscillation amplitude has a maximum value for a temperature much below 
the Thouless temperature, T <C Tthi which qualitatively distinguish these results from those ob- 
tained for a completely transparant boundary between the mesoscopic region and the reservoirs. 
Without potential barriers between the reservoirs and the mesoscopic normal region there are no 
Andreev states that can contribute to the inter-reservoir transport. With such barriers present An- 
dreev levels are well defined and long-lived. The peak of the conductance in Fig. 2 is due to resonant 
tunneling through a macroscopic number of such Andreev levels. A small amount of non- Andreev 
(normal) quasiparticle reflection at the N/S boundaries is not, as can be seen in Fig. 3, detrimental 
to the resonant tunneling effect. Rather it results in an energy shift of the position of the reso- 
nance (provided the probability for normal reflections at the two N/S boundaries are different). 
As a result the position of the maximum of the peak amplitude is shifted to a finite temperature, 

T~ Ir&'l-lr^l T Th . 



G = ^-ir-^l TT^b^x (17) 




4. Conclusions 

We have shown that resonant tunneling through Andreev levels may give rise to an excess quasi- 
particle contribution to the normal conductance of an S/N/S sample at energies much below the 
Thouless energy and with maximal amplitude when the phase difference between the superconduc- 
tors is on odd multiple of 7r. That resonant tunnelling through Andreev levels could give rise to a 
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Fig. 2. Normalized excess conductance AG(<p)/G at of the S/N/S structure shown in Fig. la as a function of the 
phase difference (f> between the two superconductors. AG(<j>) = G((j>) — G(<f> = 0), Gjv = e r (2e 2 /h)Nj_; T/T-^h = 0.1, 




0.25 



T/T Th 

0.1 0.2 0.3 0.4 0.5 



Fig. 3. Temperature dependence of the maximum normalized excess conductance AG(ij> = tt)/Gm for the same 
tunnel barrier parameters as in Fig. 2. 
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"giant" effect was proposed for the case of a ballistic normal region in Ref . || . The effect discussed 
here is not sensitive to the details of the semiclassical motion of quasiparticles inside the disordered 
normal region. A more remarkable fact is that the main conclusion about the importance of res- 
onant transmission through a macroscopically large number of Andreev levels seems to be valid 
not only for a smooth disorder potential allowing a semiclassical analysis to be carried out as done 
here. Indeed, one can consider the density of states in an S/N/S junction isolated from reservoirs 
and solve the Eilenberger-Usadel equation, which is valid for a short-range disorder potential in the 
normal region. This approach gives, according to Ref. E§], a gap in the spectrum which closes at 
4> = 7r with the density of states diverging at E = 0. 

In recent experiments |37| fll| oscillations of the conductance G((j>) have been observed, which 
are well described in the framework of the thermal effect of Nazarov and Stoof at temperatures 
around the Thouless temperature, T ~ Tj^, but are significantly difercnt from what the thermal 
effect can explain at lower temperatures. The conductance maxima were found to be 7r-shifted in 
both experiments in the range of low temperatures. The maximum oscillation amplitude AG(T) was 
observed at T w 20 mK (the Thouless temperature was w 200 mK pH). One may speculate (see 
[0) that grain boundaries and geometrical feauters of the contacts act to split the semiclassical 
quasiparticle trajectories in the same sense as the potential barriers in the model used here. If so 
one would expect the thermal effect J?], |f| and the resonant tunneling effect to co-exist. In this case 
the cross-over in the phase and temperature dependences of the conductance oscillations could be 
understood as a result of competition between the "high temperature" thermal effect and the "low 
temperature" resonant tunneling through Andreev levels. 
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Foundation under Grant No. PHY94-07194 is gratefully acknowledged. MJ is grateful for the hos- 
pitality of the Institute for Theoretical Physics, UC Santa Barbara, where part of this work was 
done. 



Appendix 1 

The set of equations (^j can be re-arranged in the following way, 

{D( £ )+e r p}|<z) = ^l<? (e) >- (18) 

Here Q(e) = 1 — U(e), 1 is a unit matrix of order 2N x 2N (N is the number of sections in the chain 
in Fig. lb), and U{e) is a 2N x 2N unitary matrix whose explicit form can be found from Eq. (||) 
by setting e r = 0. The components of the vector \a) are the amplitudes <z„ , where the subscript 
n labels the sections of the chain (n = 0, ±1, ±2, ... ± N/2); the superscript i denotes the electron 
(i — e) and hole (i = h) paths, P is an operator projecting onto the section of injection 

P=(l/2){\gV)(gW\-\gW)(g^\} . (19) 

The vector \g^) has only one non-zero component 8 n fi8 i e as has \g^), 5 n fifii,h [see EQ. (0)]. The 
energy spectrum of a quasiparticle moving along the path of Fig. lb in the absence of any coupling 
to the reservoirs is determined by the roots of the determinant of the matrix uj. In order to find 
the resonant transmission amplitude we apply resonant perturbation theory to the set of algebraic 
equations ( [l8| ) assuming the barrier transparency e r to be small (e r <C 1) and the energy E of the 
incoming qusiparticle to be close to the energy level E a . Expanding the matrix u)(E) to lowest order 
in energy, Q(E) = ui(E a ) + lj (E a ) (E — E a ) and expanding the vector |o) to lowest order in the 
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small parameter ^/e~ as 

|a)=7|e a )+V^|aW), (20) 

one gets a set of algebraic equations, 

u(E a )\a,W) = -1 {(E-E a ) Q'(E a ) + e r P}\e a )+^;\g) (21) 

where the vector \e a ) is a normalized non-trivial solution of the equation^ 

w(E a )\e a >= 0. (22) 

The constant 7 is determined by the condition that Eq. (|2l]) has nontrivial solutions. Its value 
can be found by multiplying both sides of the set of equations with (e a \, the latter vector being a 
non-trivial solution of the set of equation, 

{e a \Q(E a ) = (23) 
(the vectors (e Q | and \e a ) are normalized in such a way that (e a \e a ) = 1). As a result, one gets 

7 = V V iy ' - , (24) 

(ed uj'(E a )\e a ) (E - E a ) + e r (e a \P\e a ) 

with (e a I ) = et§ a . As we are considering the case of low normal reflection probability amplitudes, 
\rpf I <C 1, and the terms in the denominator already contain small factors (E — E a ), the amplitudes 
r^j can be neglected in these terms to lowest order. This gives as a result that (e a \Q (E a )\e a ) = i/E a 
(Eo = Hvf/to, tq is the propagation time along the path in the section of injection in Fig. lb) and 
(e a |P|e a ) = 1/2. 

One can show that in the resonant approximation the second term on the right hand side of 
Eq. (EG) can be neglected if the prameter e r -C 1. Hence the amplitude of the electron- hole transition 
is Oq = ^{g^\a). Therefore, using Eq. (f7j) and ( p^ ) we get the Breit-Wigner formula (||). 

5. Appendix 2 

As shown by Slutskin [Q, the density of states of a one- dimensional chain of the type in Fig. lb 
can be written as a Fourier series, 

1 N 

71=1 

where 

F n =J2D({ln})U eilnrmE/h - 

{lm} m 

Here {/„} labels various sets of numbers {1\12,..}n) where l n —i = 1 and l m ^ n is equal to or 1; 

(1 2) 

D({l n }) are Fourier coefficients which depend on the "hopping integrals" r N ' . Since the "Fourier 
amplitudes" do not depend on r„ one only has to average the product r n Y[ rl e ll ^ T ^ E / h while calcu- 
lating the average density of states. Using the Lorentzian distribution ( |l3| ) one finds the result 



e U n r El /h ^ 



f Here and below, there is no summation with respect to double indices. 
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The term in parenthesis in this expression is exactly the same as what appear on the right hand 
side of Eq. ( |25| ) provided the system in Fig. lb is periodic with all r„ = f. From this and the fact 
that vq(Ei) — Vq(—Ei) the result Eq. (|l4|) follows immediately. 
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